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Abstract 

Broadband inversion pulses that rotate all magnetization components 180° 
about a given fixed axis are necessary for refocusing and mixing in high- 
resolution NMR spectroscopy. The relative merits of various methodologies 
for generating pulses suitable for broadband refocusing are considered. The 
de novo design of 180° universal rotation pulses (180^^) using optimal con- 
trol can provide improved performance compared to schemes which construct 
refocusing pulses as composites of existing pulses. The advantages of broad- 
band universal rotation by optimized pulses (BURBOP) are most evident for 
pulse design that includes tolerance to RF inhomogeneity or miscalibration. 
We present new modifications of the optimal control algorithm that incor- 
porate symmetry principles and relax conservative limits on peak RF pulse 
amplitude for short time periods that pose no threat to the probe. We apply 
them to generate a set of ^^0°qurbop pulses suitable for widespread use in 
^^C spectroscopy on the majority of available probes. 
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1. Introduction 



Many NMR applications require refocusing of transverse magnetization, 
which is easily accomplished on resonance by any good inversion pulse sand- 
wiched between delays, ie, the standard A-180°-A block. For broadband 
applications, a universal rotation (UR) pulse that rotates any orientation of 
the initial magnetization 180° about a given fixed axis is required to refocus 
all transverse magnetization components. A simple hard pulse functions as 
a UR pulse only over a limited range of resonance offsets that can not be 
increased significantly due to pulse power constraints. 

Although a great deal of effort has been devoted to increasing the band- 
width of inversion pulses, most broadband inversion pulses [ll, [I, [s], 0, 0, @, 



H i 1, 0, 



11, 12, hj 



fly, 



16 



execute only a point-to-point (PP) rota- 
tion for one specific initial state, magnetization M-, — )■ — M^, and are not 
UR pulses. However, two PP inversion pulses with suitable interpulse delays 
can be used to construct a refocusing sequence [l7|, [l8|, which is effectively a 



360^^ pulse. Alternatively, a 180?^^ refocusing pulse can be constructed from 
three adiabatic inversion pulses 19| with either pulse length or bandwidth 
of the adiabatic frequency sweep in the ratio 1:2:1. More generally, we have 
shown that one can construct a UR pulse of any flip angle from a PP pulse 
of half the flip angle preceded by its time- and phase-reversed waveform j2o| . 
Thus, a 180^^ pulse can be constructed from two 90pp pulses. 

The reliance on composites of PP pulses to construct UR pulses high- 
lights the perceived difficulty of creating stand-alone UR pulses. The de 
novo design of UR pulses for NMR spectroscopy has received compara- 
tively little attention joi 21|, so it is an open question whether the com- 



posite constructions using PP pulses achieve the best possible performance. 
Yet, the demonstrated capabihties of optimal control for designing PP pulses 



22|, |23|, 124, |25|, l26|, |27|, [28 J29 
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3ll l32l . l33l l34j . |35| are equally applicable to 



the design of UR pulses [36|, |37|, l38|]. The required modifications to the basic 
optimal control algorithm are fairly straightforward 36|, l39|, |40| and maintain 
the same flexibility for incorporating tolerance to variations in experimentally 
important parameters, such as RF homogeneity or relaxation. 

In this work, we design broadband refocusing pulses by optimizing the 
propagator for the required UR transformation. The resulting ^9>Q°burbop 
pulses (broadband universal rotation by optimized pulses) are compared to 
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existing composite refocusing pulse schemes to characterize the conditions 
under which one design method might be preferrable to another. In addition, 
we introduce new optimal control strategies tailored to take advantage of 
specific opportunities available in the design of UR pulses. The culmination 
of these efforts is a set of low-power, high-performance broadband refocusing 
pulses that satisfy the power constraints of widely available probeheads and 
complex multipulse sequences. 



2. Optimal control algorithm for 180^^ pulses 

A general procedure for creating a desired unitary propagator in an arbi 



trary (closed) quantum system is given in 36|, |39|, |40|, |4l| . Time evolution 
proceeds according to a matrix exponential of the system Hamiltonian. For 
two-level systems, as in many NMR applications involving a single nonin- 
teracting spin-1/2 species, this evolution is well-known to be equivalent to 
a rotation of the 3D vector representing the state of the system about the 
effective applied field |42|. The relatively abstract general procedure for 
propagator optimization can be made considerably more transparent in this 
case. 



2.1. Flavor I (basic vanilla) 

The optimal control methodology for generating PP transformations in 
two-level systems has been described in detail previously [22|, |23|, |2J, |25|, l26|. 



271 . l28l l29l . l30l . l3ll . l32l . l33l l34j . |35| . with progressive modification to enhance 
pulse performance by incorporating various experimental constraints. In each 
iteration of the algorithm, one starts with a given initial magnetization Mq, 
applies the RF derived for the current iteration, and compares the resulting 
final state Mf with a desired target state F. This comparison, quantified 
in terms of a cost function, allows one to efficiently calculate a gradient for 
improving pulse performance in the next iteration. A simple but effective 
cost function $, which we employ in the present work, is the projection of 
the final state onto the target state, the standard inner product ( F \ Mf ) , 
given in this case by the dot product F ■ Mf- Desired RF limits can be 
enforced by clipping without degrading the performance of the optimization 



28 



The algorithm for generating UR pulses in the single-spin case is a straight- 
forward modification of the PP algorithm. A 180^^ pulse applied, for exam- 
ple, along the y-axis to magnetization M effects the transformation (M^, My, M^ 
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{—Mx, My, —Mz). This is simply three separate PP transformations of the 
initial states Mi = (1,0,0), M2 = (0,1,0), M3 = (0,0,1) to their respec- 
tive target states Fi = (-1,0,0), F2 = (0,1,0), F3 = (0,0,-1). The cost 
function comparing the final states Mkf {k = 1,2,3) at the end of an RF 
pulse to the target states is 



$ = Fi ■ Mif + F2 ■ M2f + F-i- M 



'■■if- 



(1) 



The algorithm proceeds in the standard fashion using this cost function. We 
will refer to this as algorithm A. 

This simple intuitive modification to the cost is exactly equivalent to an 
analogous procedure given in [36] for optimizing the unitary propagator, 
which can be seen as follows. The rotation operator Rp m. 3D corresponding 
to the target propagator that generates a 180° rotation about the y-axis is 
given by 




( 



F\ F2 F3 



\ 



\ 



J 



(2) 



ie, the i*^ column is the corresponding PP target F^. 

The actual rotation operator at the end of a pulse of length Tp is 



RiTp 




Ml M2 Mg 



v 



Ml/ M 



2/ 



M 



3/ 



(3) 



with the rotation transforming each column to its associated final state for 
the individual PP transformations. 

The cost is again given by the projection of the final state onto the target 
state, with the inner product 



^r={Rf\ R{Tp) ) = Tr [Rl RiTp)] 



(4) 
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where superscript T denotes the transpose, and the operator Tr returns the 
trace (sum of diagonal elements) of its argument. We then have 



RpR{Tp) 




/ ; ; ; \ 

Ml/ M2/ M3/ 



(5) 



V : : : / 

The sum over diagonal elements of this matrix product gives Eq. [[1]]. 



2.1.1. Flavor II ( symmetry principle ) 

The formalism for constructing UR pulses from PP pulses provides 
additional insight for improving ^'&^°burbop performance. The symmetry of 
the construction procedure constrains the resulting rotation axis to the plane 



defined by the desired axis and the z-axis [43|. Details of the results which 
follow are provided in the Appendix. 

For a UR rotation about the x-axis, any nonidealities in the original PP 
pulse (for example, due to resonance offset) shift the resulting rotation axis 
only in the x2-plane. Given a phase deviation of magnitude 50 from the 
desired target state and rotation error 56 compared to the desired rotation 
angle, the angle aigo^,^ that the rotation axis in the x2;-plane makes with the 
X-axis is small and bounded according to the relation 



«180, 



UR ~ 



< v/(509Opp)2+(5^^9Opp)2 (6) 



for small deviations measured in radians. 

This has the effect of reducing any phase errors in the original 90pp pulse 
used for the construction. For example, a 180° rotation about any axis in 
the a;2;-plane transforms ly to —ly with no phase error. There would be no 
phase error in the — )■ Ix transformation, either, but the amplitude of the 
final x-magnetization would decrease depending on the angle cxi^q^j^, giving 
a cos(2Q;i80y^) dependence. Similarly, the — t- — transformation would 
have a cos(ai80ujj) dependence for the magnitude of the final z-magnetization. 
Phase errors are thus the result of deviations, drj, from the ideal 180° rotation 
angle, which are bounded (see Appendix) according to 

'5r7<2V(<509Opp)2 + (5^9Opp)2 (7) 

Larger amplitude errors result from both the displacement of the rotation axis 
from the x-axis and deviations from the ideal 180° rotation angle. Phase de- 
viations in the 90pp pulse are reduced in the resulting 180^p pulse according 
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to the relation 

^<t>180un ^ (<^09Opp)' + {S9gOppf. (8) 

The reduction is quite significant, with SOgopp = 1° and 509Opp of 10°, 3°, 
and 1°, for example, giving bounds for (^^isojj^ of 1.76°, 0.17°, and 0.03°, 
respectively. 

For applications requiring high phase fidelity which can afford modest 
loss of signal intensity, we therefore incorporate the symmetry principle of 
the construction procedure into the optimal control algorithm. For RF pulse 
components and Uy digitized in time steps, the first half of the pulse is 
determined using the basic algorithm A. The second half of the pulse is then 
constructed using the time-and phase-reversed components of the first half. 
Phase zero for Ux leaves it unaffected, while Uy is inverted to give 

^ _„^/2+l-i (9) 

for i = 1, 2, 3, ... , N/2. We refer to this algorithm incorporating the symme- 
try of the construction principle as A5. 

2.1.2. Flavor III (time- dependent RF limit) 

Peak RF amplitude must remain below probe limits (e.g., available for 
^^C spectroscopy), but larger RF amplitudes result in improved broadband 
performance. For sufficiently short time periods, we note that probe RF 
limits can be higher than conservative limits that protect the probe from 
arcing under any conditions. Enforcing a lower probe limit for an entire pulse 
duration can sacrifice performance unnecessarily. We therefore introduce a 
time-dependent RF field limit to allow increased RF amplitude for short time 
intervals and achieve improved performance. Empirically, we find that low 
RF limits force algorithm A to request higher RF amplitude in the middle 
of the pulse for improved performance. We therefore allow a higher RF limit 
for a short time during the middle of the pulse. We refer to this algorithm 
as At, or, if it is also combined with the symmetry principle, as Ast- 



3. BURBOP compared to refocusing with PP pulses 

Performance of the new ISO^^/^^op pulses designed using optimal control 
are compared with previous methods for generating refocusing pulses, start- 
ing with those that can be constructed from 90pp pulses. We then consider 
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the refocusing performance of two 180pp pulses. In both these comparisons, 
RF amphtudes are hmited to a relatively conservative peak value of 11 kHz 
(22.7 /is hard 90° pulse) for widespread use in ^^C spectroscopy. Finally, 
we present a set of ISO^^yjij^Qp pulses that utilize the higher RF power lim- 
its allowed for short time intervals to increase the maximum RF amplitude 
from 11 kHz to 15 kHz (16.7 hard pulse) in the middle of the pulse. The 
performance of these pulses is compared with the composite adiabatic pulse 



scheme [19| implemented as a smoothed Chirp pulse [44, |45| in standard 



Bruker software 



3.1. Algorithms A and As 
3.1.1. Two 90°pp pulses 

A previously published 90pp pulse with exceptional performance [s^] was 
used to construct a 180^p pulse according to the procedure in |2^ . This 1 ms 
constant amplitude phase-modulated excitation pulse transforms greater than 
99% of initial z-magnetization to the x-axis over a resonance offset range of 
50 kHz for RF amplitude anywhere in the range 10-20 kHz. For all but 
the lowest RF amplitudes, the transformation is greater than 99.5%. Phase 
deviations of the excited magnetization from the x-axis are less than 2°-3° 
over almost the entire optimization window, with minor distortions in the 
6°-9° range at the lowest RF values. 

The upper panel of Fig.[T] depicts the performance of the constructed pulse 
composed of two 1 ms 90pp pulses. As expected from the discussion in §2.1.H 
the signal intensity of the transformed magnetization is decreased slightly 
compared to the original PP pulse, but the phase relative to the target axis is 
improved considerably. The construction procedure rotates phase deviations 
relative to the target axis out of the transverse plane, reducing spectral phase 
errors that are the primary interest. Increased phase dispersion is allowed in 
a plane orthogonal to the x-y plane, resulting in reduced projected transverse 
magnetization. 

In comparison, the basic algorithm A, which optimizes the UR propagator, 
results in a pulse with performance shown in Fig.llb that is more consistent 
with the performance characteristics of the 90pp pulse used for the UR con- 
struction in Fig. [1^. Thus, algorithm A maximizes signal amplitude, with 
phase fidelity secondary, while the construction procedure maximizes phase 
fidelity, with signal amphtude secondary. 

When the optimization includes a range of RF inhomogeneity/miscalibration, 
algorithm A can produce shorter pulses than the construction principle for a 
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comparable amplitude performance, as shown in Fig.Ht, at the cost of poorer 
phase performance. We find empirically that the reduction in pulse length is 
50%, 30%, and 15% for RF miscalibration ranges of ±33%, ±10%, and ±5%, 
respectively. In optimizations that do not include tolerance to RF inhomo- 
geneity, there is no reduction in pulse length compared to the construction 
procedure. 

Finally, algorithm As utilizes the construction procedure symmetry in 
optimizing the UR propagator, producing the pulse performance shown in 
Fig.[T]i. Algorithm provides more ideal phase performance and therefore 
can provide an advantage over the construction procedure when tolerance 
to RF inhomogeneity is included in the optimization. The choice between 
algorithms A and As depends on the application, whether maximizing signal 
amplitude is more important than ideal phase performance. 

3.1.2. Two 180°pp pulses 



Two 180pp pulses can be used to refocus transverse magnetization [17|, 



18| by modifying the standard spin-echo sequence. For example, a 90pp 
excitation pulse of length Tp that produces a linear phase dispersion AuRTp 
(0 < -R < 1) as a function of offset Aco could be followed by r - 180pp 
- (r ± RTp) - 180pp. The first r-delay already includes a phase evolution 
equivalent to time RTp, which is included in the next delay of the standard 
spin-echo. The second inversion pulse then compensates the phase errors of 



the first [14[, which is unnecessary if a single 180^p is employed. A shorter 
sequence for this example is (r — RTp) - 180pp - r - 180pp. The resulting 
modified spin-echo procedures rotate magnetization 360° and are not UR 
inversions, but this is not an issue for refocusing. 



BIP M 

and BIBOP [15|, |16[ pulses are optimized to provide exceptional 
performance as 180pp pulses, which translates to similar performance when 
incorporated into the above refocusing scheme. There is little room for im- 
proving the refocusing performance of two 180pp pulses. A similar conclusion 



was reached in the context of relatively high bandwidth selective pulses |46 
On the other hand, a single 180^^ pulse can be simpler to incorporate into 
complex pulse sequences with respect to adjusting the timings and synchro- 
nization among various pulses. Conventional hard 180° pulses can be easily 
replaced in an existing sequence by 180p^^pQp pulses without further change 
of the pulse sequence. By contrast, incorporating two 180pp pulses can re- 
quire the adjustment of phases and a correspondingly detailed understanding 
of the pulse sequence. 



8 



We therefore find at most a modest advantage in using algorithms A or A5 
to generate a 180^j:j pulse compared to the overlapping spin-echo sequence 
using optimized 180pp pulses. However, incorporating a time- dependent RF 
limit into the optimal control algorithm does provide a distinct advantage for 
generating a single 180^^ pulse compared to existing pulses, as illustrated in 
what follows. 

3.2. Algorithms At and As^t 

Broadband refocusing bandwidths of ~ 50 kHz, sufficient for high-field ^^C 
spectroscopy, are readily achieved using any of the pulse schemes discussed 
so far. The peak RF power required for the pulses is well-within hard-pulse 
power limits for modern high resolution probes. However, in multipulse se- 
quences, repeated application of what might be deemed a modest power level 
for a single pulse can be a problem if the total energy delivered to the sample 
(integrated power) is too high. There are also limits on the total energy 
that can safely be delivered to a given probe. For these reasons, the most 
general and widespread applications impose peak power levels that are more 
conservative than what might be necessary for broadband refocusing using a 
typical probe. We therefore incorporate a time- dependent RF limit into the 
optimal control algorithms to keep peak power low for most of the pulse, but 
allow short increases in this limit where it can have the most benefit. We 
utilize algorithms At and As^t (defined in §2.1.2p to investigate the possibil- 
ity of generic broadband refocusing pulses suitable for use with any standard 
probehead in any pulse sequence. 

3.2.1. Three adiabatic inversion pulses 

The best broadband refocusing performance available in the standard 
Bruker pulse library satisfying the required conservative pulse power limits is 
obtained using the pulse designated Chirp80 |45|]. This pulse is constructed 
from three adiabatic inversion pulses with pulse lengths in the ratio 1:2:1 



19[. It utilizes for its shortest element a 500 ms smoothed chirp pulse [44 
with 80 kHz sweep. The first 20% of the pulse rises smoothly to a maximum 
constant RF amplitude of 11.26 kHz according to a sine function before de- 
creasing in the same fashion to zero during the final 20% of the pulse. The 
final pulse is thus 2 ms long. 

Maintaining this pulse length and mindful of the given conservative peak 
RF amplitude, we designed the set of four pulses listed in Table 1. For 
most of the pulse, the nominal RF amplitude is a constant 10 or 11 kHz. A 
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maximum RF amplitude of 15 kHz is applied for 60 /is in the middle of the 
pulse, as illustrated in Fig.HJ This short increase in pulse amplitude provides 
significant improvement in pulse performance compared to ChirpSO. 

The amplitude profile shown in Fig. |2] is reminiscent of the hyperbolic 
secant pulse joj , which maintains a low amplitude for most of the pulse with a 
peak in the middle. All four pulses show excellent performance over the listed 
ranges in offset and RF field inhomogeneity/miscalibration. Performance is 
comparable to the performance shown in Fig. [1] for higher power, constant 
amplitude pulses with nominal peak RF of 15 kHz. Pulse 4 provides the 
most relevant comparison, since it has a similar range of tolerance to RF 
inhomogeneity. As expected from the earlier results for algorithms A and 
As, the best amplitude performance is obtained by algorithm and the 
best phase performance by algorithm As^t 

Figur^ compares theoretical performance of pulses 1 and 4 from Table 1 to 
the ChirpSO pulse. The new pulses significantly improve phase performance 
over the targeted range of offsets and RF inhomgeneity/miscalibration. Ad- 
ditional quantitative comparison between pulse 4 and ChirpSO are provided 
in Figs JH and which also show the excellent agreement between simulations 
and experimental pulse performance. Improvements in lineshape and phase 
that are possible using the new pulses are shown in Fig. [61 



Pulse 


Algorithm 


R-Pnominal 


RFmax 


Inhomogeneity 


Transformat: 


ion Error 






(kHz) 


(kHz) 


optimization 


(amplitude) 


(phase) 


1 


At 


11 


15 


±10% 


< 0.2% 


< 2° 


2 




11 


15 


±15% 


< 0.9% 


< 0.3° 


3 




10 


15 


±15% 


< 1.5% 


< 0.5° 


4 




11 


15 


±25% 


< 2% 


< 0.8° 



Table 1: Four pulses optimized to execute a 180° universal rotation about the y-axis over 
resonance offsets of 50 kHz and RF field inhomogeneity listed in column 5. All pulses are 
2 ms long. For all but 60^s, pulse RF is constant at the value RFnominai, increasing to 
RFmax in the middle of the pulse (see Fig.[2|). Amplitude errors in the transformation 
are the maximum deviation from the target magnetization over the optimized ranges of 
resonance offset and RF inhomogeneity, expressed as a percent. Similarly, phase errors 
represent the maximum deviation from the target phase. Performance of pulses 1 and 4 
are shown in more detail in Fig. [3] 
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4. Experiment 

All experiments were implemented on a Bruker 750 MHz Avance 111 spec- 
trometer equipped with SGU units for RF control and linearized amplifiers, 
utilizing a triple-resonance PATXl probehead and gradients along the 2;-axis. 
Measurements are the residual HDO signal in a sample of 99.96% D2O doped 
with CUSO4 to a Ti relaxation time of 100 ms at 298° K. Signals were obtained 
at offsets between —25 kHz to 25 kHz in steps of 500 Hz. To demonstrate the 
tolerance of the pulses to RF inhomogeneity/miscalibration, the experiments 
were repeated with RF amplitude incremented by ±15% and ±25% relative 
to the nominal maximum RF amplitude for each pulse (15 kHz for pulse 4 
of Table 1, 11.26 kHz for ChirpSO). To reduce the effects of RF field inho- 
mogeneity within the coil itself, approximately 40 /zl of the sample solution 
was placed in a 5 mm Shigemi limited volume tube. 



5. Conclusion 

We have presented three optimal control algorithms for the de novo de- 
sign of universal rotation pulses, applied specifically to inversion. The most 
noteworthy innovations for NMR spectroscopy are inclusion of the construc- 



tion principle discovered in [20[ and a time- dependent or "floating" limit 
on the peak RF allowed during the pulse. The new algorithms result in 
improved performance compared to existing UR pulses constructed as com- 
posites of point-to-point pulses. The methodology is very general, and fur- 
ther improvements in the design of robust universal rotation pulses can 
be anticipated. The new pulses are implemented in the Bruker pulse li- 
brary and will also be made available for downloading from the website 



http: / / www.org.chemie.tu-muenchen.de/glaser/Downloadsiitml , 
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6. Appendix 

The results of section 12.1.11 were derived as follows. 

6.1. Euler angles 

The rotation of a vector can be represented as the composition of three el- 
ementary (Euler) rotations about the fixed axes used to represent the vector: 
rotation by angle ip about the z-axis followed by rotation 9 about the x-axis 
and rotation ip about the z-axis. Using the notation cp = cos/3,S/3 = sin/3, 
and Rk{(3) for a rotation by angle /3 about axis k, we have 

R^ii^) = 
RM = 

R,{^) = \ ce -se \ (10) 







-s^ 














1 


1 











Ce 


-S0 





Se 


Ce 


1 











Ce 


-Se 





Se 


Ce 



givmg 



R = R,{v)RMRzW 




Cipc^ ceSnjjS^ CipSip Cr^ceSfp 
cec^Sip ~\~ CfpSip Cfpcec^ s^s^ 
SipSe c^se 

Rii Ri2 Ri3 \ 

R21 R22 R23 (11) 

-R31 R32 R33 / 

The Euler angles describing the rotation are thus 

9 = COS-l[i?33] 

if = tan-i[/?i3/(-/223)] 

^ = tan-i[i?3i/i?32] (12) 

Consider the rotation matrix W that produces the 90pp transformation 
Iz — > —ly For initial Iz, the first z-rotation has no effect. Angle 9 therefore 
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gives a direct comparison with the desired rotation angle of 90°, while gives 
the phase deviation from the target state —ly. W also rotates initial states 

and ly to final states that can be arbitrary in the design of the original 
PP transformation. This provides a great deal of flexibility in designing PP 
transformations compared to UR transformations, which must rotate each 
initial state to a specific final state. For an ideal ^^°uji rotation about the 
X-axis, Tp would be zero. It thus provides the orientations of the rotated 
vectors WIx and Wly in the plane orthogonal to the final state WIz- 

The pulse that generates a 180^j:j transformation from this particular 90pp 
pulse is constructed from the phase-inverted pulse, denoted by 90pp, followed 



by the time-reversed (tr) pulse 90pp 20|. At a resonance offset —v, we have 
the corresponding rotation operators 

W''\-v) = Rz{l80)W-\u) Rz{180) 

Wu W22 I (13) 

and 

W{-iy) = Rx{180)W{iy) Rx{180) 

-W21 W22 W2Z I , (14) 

-W^31 W^32 W^33 

where /^^(ISO) is zero except for (—1,-1,1) along the diagonal, -Ra;(180) is 
similar, but with (1, —1, —1) along the diagonal, and is the transpose 
oiW. _ _ 

The 180^j:j rotation matrix W{—v) = W^"^ {—u)W {—u) therefore differs 
from the identity matrix to the extent terras^ in the product for each matrix 
element have different signs. For example, Wu is equal to W^i — W21 + W^i 
rather than W^f^ + W21 + W^^ = 1. Adding and subtracting W21 in the 
expression for Wu gives 1—214^2^ Preceding similarly, all the matrix elements 
can be written in terms of their difference from zero (off-diagonal) or ±1 
(diagonal) to obtain 

1 - 2Wi^ 2W21W22 2W2xW2z 

W = \ -2Vr2lVr22 -1 + 2iy|2 21^23^^22 | (15) 

21^21^^23 -2W^22W^23 -1 + 2{W^s + W^) 
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for the 180^^ rotation operator W at resonance offset —u in terms of elements 
Wij of the 90pp rotation operator at offset u. Differences in performance 
at offsets ±z/ are small to the extent the PP optimization is successful in 
generating uniform performance over a symmetric range of resonance offsets. 

The Euler angles for W obtained using Eq. [ [T^ ] give (p = if). This result 
provides geometric insight into the performance of W . Rotation about the 
2;-axis followed by 180° rotation about x then rotation by the original amount 
about z will return x to x and send y and z to —y and —z, respectively. The 
deviation from the ideal goal of a 180° rotation determines the errors in the 
180^^ transformation. For initial I^-, the phase angle for the rotated state 



WIx relative to the target state is tan0 = W21/W11. Using Eq. [ [TT] ] and 
adapting our earlier notation to = cos/3, etc., gives 



tani 



c^s^ (1 + Ce) 



tan ■ 



Ce tan^ ip 



(16) 



The extrema for tan0 as a function of ijj occur for dltan (p) / dip = 0. There 
are two solutions, one of which is an inflection point and the other producing 
a maximum for tanip = —Cn^^'^, so that 



tan( 



l + ce 



2V^e 



(17) 



According to Eq. [ [T2] ] and Eq. [[T5]]. 



cos6' 



33 



-1 + 2(1^1^3 + 



33 J 



-1 + 2(^^4 + 0^) 



-1 



(18) 



where we have substituted 9 = n/2 + 69. 

Both 69 and are small for a good 90pp pulse. Expanding cse ^ ^ — 
{69)'^/2, Sip ~ 99, s^e ~ 69, performing the multiphcations in Eq. [[TSl]. and 
keeping terms to second order gives, with 9 = + 69, 



cos 9 



-1 + 2(69^ + ^^ 
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COS (5^ ^ l-2{5e^ + ^^) 

56 ^ 2^W+^. (19) 

Thus, the deviation 59 from the ideal 180^j:j rotation angle is slightly more 
than twice the deviation 59 from the ideal 90pp rotation angle. 

Using Eq. [[19]] for cos^, the small-angle approximation for tan0, and 
applying (1 — x)^^"^ ~ 1 —x/2 for small x to tan?/; = l/(— cos^ )"'^/^ gives 

t ^ < + 

^ (5^2 ^¥^2 (20) 

Similarly, for initial ly, the phase relative to the target — y-axis is tan0 = 
W12/W22, which results in tanip being replaced by cotip in Eq. [US]]. The 
solution for tan0max then occurs for cot-?/^ = (— cos6')~^/^, giving the same 
bound for 0max as in Eq. 



6.2. Angle/axis parameterization 

Alternatively, we can describe a rotation in terms of rotation angle t] about 
an axis defined by unit vector n to obtain 

_i Tr i? - 1 

rj = cos 



Uy 



2 

R32 — R23 

2 sin 77 
-R13 — R31 

2 sin rj 



R2I — R12 

i^z = -1^. 21 

2 smrj. 

Applying this to W in Eq. [[15]] shows = 0, as expected from symmetry 
arguments noted in section 12.1.11 The rotation axis^ makes an angle orwith 
respect to the x-axis given by tana = Uz/ux = (W21 — Wi2)/(W32 — W23). 
Following arguments similar to those leading to Eq. [[16]] gives, for small a, 

_ 2c^s^ (1 + ce) 

2c^se 
< {l + ce)/Se 

^ ^fWT^. (22) 



a ~ tan a 
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Here, we have written sq = (1 — c^)^/^, substituted from Eq. [ IT^ ] keeping 
terms to second order, and used the maximum value of one for s^. 

We obtain similarly from Eq. [ flVi ] cos 77 = — 1 + 214^2^2 ~ ^22 after rear- 
ranging terms as in the discussion preceding Eq. [[T5]]. Substituting as above 
gives W22 = c^ce — s^. Writing = 1 — and employing the maximum 
value of one for gives 

cos?7 < cos 6* (23) 

and the deviation of each angle from the ideal rotation of 180° can be rela- 
tively quantified as 

6r] < 56, (24) 
with 66 given in terms of the 90pp angles in Eq. [[T9]]. 
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Figure 1: Theoretical performance of four 180^^ pulses for inversion of magnetization 
about the y-axis, designed using different algorithms. The first two columns show the 
transfer efficiency for the labeled transformations relative to ideal or complete transfer. 
The last column displays phase deviation in degrees relative to the labeled target state. 
The nominal peak RF amplitude for all th^jjulses is i3j = 15 kHz, optimized to perform 
over a resonance offset range of 50 kHz and variation in RF homogeneity /calibration of 
±33%. All pulses are constant amplitude with the exception of pulse b), which deviates 
from the maximum for less than 10% of the pulse, a) constructed from the 1 ms 90pp 
pulse of Ref. [s^] preceded by its time- and phase-reversed waveform [i^, pulse length 
Tp = 2 ms. b) algorithm A, Tp = 2 ms. c) algorithm A, Tj, = 1 ms. d) algorithm As, 
which incorporates the symmetry principle used in a), Tp = 2 ms. 
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Figure 2: The amplitude and phase of 180^^^. pulse 4 from Table 1 obtained using algorithm 
As,T- A conservative limit of 11 kHz for the peak RF applied for 2 ms is relaxed to allow a 
safe peak of 15 kHz for 60/us. The pulse is amplitude-symmetric and phase-antisymmetric 
in time, incorporating the symmetry of the construction procedure from Ref. [20| . 
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Figure 3: Similar to Fig.[Tl but the maximum RF amplitude is allowed to float for a 
sufficiently short time period, giving the reduced-power pulse in Fig. [2] All three pulses 
are designed to operate over a resonance offset of 50 kHz and a range of variation in RF 
homogeneity/calibration relative to the ideal B^, as given in Table 1. a) pulse 1 of Table 
1, RF tolerance ±10%. b) pulse 4 of Table 1, RF tolerance ±25%. c) composite adiabatic 
refocusing [l^ using pulse ChirpSO from the Bruker pulse library (isj . 
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Figure 4: Further quantitative detail for the —J' —M^ transformation from Fig.[3}D for 
pulse 4 of Table 1 (black) and Fig.[3t for ChirpSO (red), plotted for RF scalings of ±15% 
and ±25% relative to the nominal maximum RF amplitude at 0% (15 kHz, pulse 4 and 
11.26 kHz, ChirpSO). Theoretical values for the inversion profile are plotted on the left as 
a function of resonance offset, with phase deviation ip relative to the target —Mx plotted 
on the right. Adiabatic ChirpSO produces significant phase errors within the bandwidth 
for all RF scalings, in contrast to the almost ideal performance of the optimal control 
^^%URBOP pulse. 
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Figure 5: Experimental measurements of the inversion profile (left) and phase deviation 
(right) corresponding to the simulations in Fig.lU showing excellent agreement between 
the experimental and theoretical performance of the pulses. 



25 



Adiabatic 



Optimal control 




RF scaling -25% -15% 0% 15% 25% 

Figure 6: Experimental lineshapes comparing optimal control 180^^ pulse 4 of Table 1 
and adiabatic ChirpSO at a resonance offset of —25 kHz for RF scalings of ±15% and 
±25% relative to the nominal maximum RF amplitude of 15 kHz at 0%. Simulations (not 
shown) are in excellent agreement with the experiments, as already shown in Figures|4] 
andO ChirpSO produces significant experimental phase errors of 20°, 15°, 9°, 11°, and 
5° (reading left-to-right across the RF scalings in the figure) in contrast to the experimental 
performance of 0.7°, 0.7°, 0.9°, 2.6°, 2.5° for the optimal control ISO^urbop P^lse. 
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